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Abstract 

We  consider  heat  balance  in  the  catalyst  layer  of  a  low-temperature  fuel  cell.  Based  on  the  heat  transport  equation  in  the  layer,  the  exact  boundary 
condition  for  the  problem  of  heat  transport  in  the  cell  is  obtained.  The  limits  of  validity  of  the  resulting  expression  are  discussed. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  rates  of  virtually  all  transport  and  kinetic  processes  in 
low-temperature  fuel  cells  (PEFCs  and  DMFCs)  exponentially 
depend  on  temperature.  Local  overheat  by  several  Kelvin  de¬ 
grees  may  change  the  rate  of  electrochemical  reactions,  shift 
liquid-vapor  equilibrium  and  enhance  diffusion  transport  of  re¬ 
actants  and  products  [1],  Clearly,  heat  balance  is  an  essential 
issue  in  fuel  cell  operation.  Proper  heat  management  is  even 
more  difficult  in  fuel  cell  stacks,  where  dense  packing  of  cells 
retards  heat  removal. 

Cell  warming  is  due  to  heat  generated  in  the  electrochemical 
reactions  and  due  to  the  electric  power  dissipated  by  currents. 
Fuel  cell  “sandwich”  supports  two  types  of  currents:  electron 
current  transported  through  the  carbon  threads  in  the  catalyst 
and  backing  layers  and  proton  current  in  the  bulk  membrane 
and  in  the  membrane  phase  dispersed  in  the  catalyst  layer  (CL). 

In  the  CL,  the  electron  conductivity  of  the  carbon  phase  is 
several  orders  of  magnitude  larger,  than  the  proton  conductivity 
of  the  electrolyte  phase.  Since  maximal  electron  and  proton  cur¬ 
rent  densities  in  the  CL  are  equal,  the  vast  majority  of  electric 
heat  there  is  released  in  the  electrolyte  phase. 

Furthermore,  typical  proton  conductivity  of  the  CL  is  at  least 
several  times  lower  than  the  conductivity  of  bulk  membrane  [2,3] 
and  thus  the  dissipation  of  electric  power  in  the  CL  exceeds  that 
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dissipation  in  the  bulk  membrane.  Taking  into  account  that  the 
reaction  heat  is  released  also  in  the  CLs  we  conclude,  that  these 
layers  are  the  major  source  of  heat  in  the  fuel  cell.1 

Numerical  models  usually  ignore  the  details  of  heat  trans¬ 
port  in  the  CLs  and  assume  that  these  layers  are  infinitely  thin 
interfaces  generating  heat  [5-7,1].  This  approach  requires  the 
expression  for  the  heat  flux  from  the  CL,  which  is  used  as  the 
boundary  condition  for  the  problem  of  heat  transport  in  the  other 
parts  of  the  cell  (backing  layers,  flow  fields,  etc.). 

In  this  work  we  consider  the  sources  of  heat  in  the  CL.  We 
show  how  heat  balance  equation  in  the  layer  can  be  reduced  to 
the  boundary  condition  for  the  problem  of  heat  transport  in  the 
fuel  cell.  This  procedure  gives  exact  expressions  and  makes  it 
possible  to  establish  their  limits  of  validity. 

2.  Heat  balance  in  the  catalyst  layer 

Consider  for  definiteness  the  cathode  catalyst  layer  (CCL) 
of  a  hydrogen  or  direct  methanol  fuel  cell  (Fig.  1).  The  “ther¬ 
modynamic”  heat  (J  mol  1 )  generated  in  the  oxygen  reduction 
reaction  (ORR)  is  given  by 

0orr  =  TAS  (1) 


1  Electron  current  in  the  backing  layers  tends  to  concentrate  near  the  edges 
of  the  current  collector  ribs  [4].  This  concentration  may  lead  to  local  overheat. 
However,  these  effects  are  out  of  the  scope  of  the  present  work. 
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Nomenclature 

b  Tafel  slope  (V) 

cp  specific  heat  ( J  g- 1  K- 1 ) 

C  molar  concentration  (mol  m-3) 

Cref  reference  molar  concentration  (mol  m-3) 

F  Faraday  constant  (9.6495  x  104Cmol-1) 

i*  exchange  current  density  per  unit  volume  (Am-3) 

je  local  electron  current  density  in  the  catalyst  layer 
(Am-2) 

/p  local  proton  current  density  in  the  catalyst  layer 
(Am-2) 

jo  current  density  in  the  cell  (A  m-2) 
y*  characteristic  current  density  (Am-2) 

/  thickness  of  the  catalyst  layer  (m) 

n  number  of  electrons  transferred  in  the  electro¬ 

chemical  reaction  (n  —  4  for  ORR  and  n  —  6  for 
methanol  oxidation) 

qj  rate  of  heat  production  due  to  currents  (W  m-3) 

qs  rate  of  thermodynamic  heat  production  due  to  en¬ 

tropy  change  (W  m-3) 

qtl  rate  of  heat  production  in  the  electrochemical  re¬ 
action  due  to  proton  transfer  through  the  potential 
barrier  (W  m-3) 

Q  total  thermodynamic  heat  produced  in  the  reaction 

(J  mol-1) 

R  rate  of  the  electrochemical  reaction  (A  m-3) 

A  S  entropy  change  in  the  electrochemical  reaction 
(J  mol-1  K-1) 

T  temperature  (K) 

x  coordinate  across  the  catalyst  layer  (m) 

~  marks  dimensionless  variables 

Greek  symbols 

a  dimensionless  parameter 

x]  polarization  voltage  of  the  CL  (overpotential,  V) 

Aeff  thermal  conductivity  of  the  catalyst  layer 

(Wm-1  K-1) 

ere  electron  conductivity  of  the  carbon  phase 

(f2-1m-1) 

<jp  proton  conductivity  of  the  electrolyte  phase 
(f2-1m-1) 

cpt  carbon  phase  potential  in  the  catalyst  layer  (V) 

( pp  electrolyte  phase  potential  in  the  catalyst  layer  (V ) 

Subscripts 

e  carbon  phase  (electron-conducting) 

ORR  oxygen  reduction  reaction 

p  electrolyte  phase  (proton-conducting) 

0  at  the  membrane  surface  (at  x  =  0) 

*  characteristic  value 

Superscripts 
w  water 

ox  oxygen 


X  j  l 


Fig.  1.  Sketch  of  the  cathode  catalyst  layer  and  the  shapes  of  electron  je  and 
proton  jp  current  densities.  Note  that  yp(0)  =  jo  is  the  total  current  density  in 
the  cell. 


where  AS  (J  mol-1  K-1)  is  the  change  in  entropy  in  ORR.  Gen¬ 
erally,  AS  depends  on  pressure  and  temperature  [13].  We  will 
assume  that  pressure  variation  due  to  the  reaction  is  negligible.2 
At  constant  pressure  the  change  in  entropy  due  to  the  ORR  is 

[13]: 


A  S(T)  =  A  So  + 


-  ASq  +  (c^  —  c°x)  In 


(2) 


where  A  So  =  A.ST/'o)  is  the  change  in  entropy  at  the  standard 
temperature  To,  Cp  and  c°x  are  specific  heats  of  water  and  oxy¬ 
gen,  respectively.  Note  that  due  to  small  mass,  the  contribution 
of  protons  and  electrons  to  the  second  term  in  (2)  is  negligible. 

Simple  estimate  shows  that  in  the  typical  for  fuel  cell  op¬ 
eration  temperature  range  (300-380  K)  the  second  term  on  the 
right  side  of  (2)  is  much  smaller,  than  the  first  one.  Physically, 
the  change  in  entropy  due  to  creation  of  new  molecules  (water) 
is  much  larger,  than  the  change  in  entropy  due  to  species  heating. 
In  the  following  we  thus  will  take  AS  =  A  So  =  const. 

The  local  rate  of  thermodynamic  heat  production  qs  (W  m-3) 
in  the  catalyst  layer  is 

TAS 

qs  =  —rprR  orr  (3) 

4  F 

where  Rorr  is  the  local  rate  of  electrochemical  reaction  (the 
proton  charge  consumed  in  the  reaction  per  second  per  unit  vol¬ 
ume,  Am-3). 

/(orr  is  the  rate  of  conversion  of  proton  current  into  electron 
current.  Generally,  Rorr  is  a  function  of  x  [8,9].  The  equation, 
which  describes  the  decay  of  proton  current  along  x  (Fig.  1)  is 

—  -Rorr  (4) 

dx 


2  The  electrochemical  reactions  in  cells  considered  occur  on  the  surface  of  cat¬ 
alyst  particles  immersed  in  Nation®.  To  our  knowledge  the  variation  of  pressure 
in  polymer  electrolyte  due  to  the  reactions  has  not  been  studied  in  literature. 
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Integrating  (4)  from  x  =  0  to  x  —  I  and  taking  into  account  that 
jp(J)  =  0,  we  get 

Rorr  dx  =  jo  (5) 

where  /  is  the  CL  thickness,  jo  =  yp(0)  is  the  local  current  in  the 
cell.  Note  that  Eq.  (5)  expresses  the  conservation  of  charge  and 
thus  it  is  valid  regardless  of  the  actual  relation  for  Rorr-3 
ORR  occurs  at  the  catalyst  surface,  which  has  potential  (pe 
(the  potential  of  electron-conducting  phase  in  the  CL).  Protons, 
however,  move  in  the  polymer  electrolyte  phase,  which  has  po¬ 
tential  (pp .  To  reach  the  catalyst  surface  the  proton  has  to  be 
transported  through  the  double  layer  at  the  catalyst/electrolyte 
interface,  i.e.,  through  the  potential  jump  —  (pp  —  cpe  (overpo¬ 
tential).  The  power  dissipated  in  this  process  is 


qn  —  q  Rorr 


(6) 


Another  source  of  heat  in  the  CL  is  electric  power  dissipated 
due  to  currents.  Since  CL  supports  electron  and  proton  currents, 
the  total  dissipated  power  is 


.  3<A 

qj  =  -Je^x 


+  Jp 


dq)p 

dx 


(7) 


where  je  and  /p  are  electron  and  proton  current  densities,  re¬ 
spectively.  Here,  we  assume  that  heat  transfer  between  the  two 
phases  is  infinitely  fast. 

Taking  into  account  Ohm’s  law  je  =  —aed(pe/dx,  jp  = 
<7pd(pp/dx  we  may  rewrite  (7)  as 


9j  = 


(8) 


Here,  <rc,  ap  are  electron  conductivity  of  the  carbon  phase  and 
proton  conductivity  of  the  electrolyte  phase,  respectively. 

Velocity  of  gas  flow  in  the  CL  is  very  small.  This  allows  us 
to  neglect  convective  term  in  the  heat  transport  equation.  Heat 
balance  in  the  CL  is,  therefore,  given  by 


3 

dx 


=  qs  +  q,,  +  qj 


or,  taking  into  account  (3),  (6)  and  (7): 


3 

dx 


J?ORR  + 


Je 


(9) 


Here,  Leff  is  the  thermal  conductivity  of  the  CL  (W  m_  1  K“ 1 ). 
This  is  the  general  equation,  which  determines  heat  transport  in 
the  cathode  CL  of  a  hydrogen  or  direct  methanol  fuel  cell. 


3  Usually,  for  Rorr  the  Tafel  equation  is  used: 

*orr=,*(^)  expG) 

Here,  i.  is  the  volumetric  exchange  current  density  (A  ni  '3),  Cox  and  Cref  local 
oxygen  concentration  and  the  reference  oxygen  concentration,  respectively,  i]  = 
ipp  —  ipe  overpotential,  <pc  and  i potentials  of  electron  and  electrolyte  phase, 
respectively,  b  is  Tafel  slope.  However,  for  the  present  analysis  the  explicit 
expression  for  Rorr  is  not  required. 


Eq.  (9)  does  not  take  into  account  the  heat  due  to  phase  trans¬ 
formation  of  water  in  the  CL.  This  process  can,  in  principle, 
be  accounted  for  as  described  in  Refs.  [10,11],  though  evapo¬ 
ration/condensation  requires  to  consider  mass  transfer  between 
the  liquid  and  vapor  phases  [12]. 

The  derivation  of  Eq.  (9)  employs  only  the  general  laws  of 
heat  and  charge  conservation.  The  only  detail  in  (9)  specific  to 
the  ORR  is  stoichiometry  coefficient  4 F.  Clearly,  Eq.  (9)  can  be 
generalized  to  describe  the  heat  balance  in  any  CL  as  follows: 


R  /2 
R+^  +  ^ 

Op 


oe 


(10) 


Here,  R  is  the  local  rate  of  the  electrochemical  reaction,  A  S  the 
respective  entropy  change  and  n  is  the  number  of  electrons  trans¬ 
ferred  in  the  reaction  (n  =  4  for  ORR  and  n  =  6  for  methanol 
oxidation). 


3.  Reduction  to  the  boundary  condition 


Calculations  show  that  in  a  wide  range  of  parameters  the 
variation  of  temperature  across  the  CL  is  vanishingly  small  (see 
Appendix  A).4  Thus,  in  many  situations  an  attractive  option 
would  be  to  replace  the  CL  with  a  thin  interface  generating  heat. 
Formally  this  means  that  in  the  “external”  problems  (e.g.  in  the 
problems  of  heat  transport  in  the  backing  layers,  flow  fields,  etc.) 
the  active  layer  is  replaced  with  the  boundary  condition,  which 
represents  the  total  heat  flux  emitted  from  the  CL.  This  approach 
is  widely  used  in  numerical  simulations;  however,  the  relations 
for  heat  flux  used  in  many  works  have  no  proper  justification. 
Moreover,  the  limits  of  validity  of  these  relations  are  usually  not 
clear. 

The  exact  boundary  condition  can  be  obtained  from  Eq.  (10) 
using  the  following  assumptions.  If  the  concentration  of  feed 
molecules  does  not  vary  significantly  across  the  CL  and  if  cell 
current: 


JO  <sC  J* 


(11) 


where 
J* 


2  Gpb 

l 


(12) 


then  the  electron  and  proton  current  densities  vary  linearly  with 
the  distance  across  the  CL  (Fig.  1): 


Jp  =  Jo 


Je  —  JO 


(13) 

(14) 


Here,  b  is  the  Tafel  slope  and  crp  is  the  proton  conductivity  of  the 
CL  (for  the  detailed  discussion  please  see  Ref.  [8]).  Furthermore, 
if  (11)  is  fulfilled,  q  does  not  vary  significantly  across  the  CL 
and  we,  therefore,  may  put  in  (10)  q  ~  const. 

Temperature  variation  across  the  CCL  is  small  (see  Appendix 
A);  thus  on  the  right  side  of  Eq.  (10)  we  put  T  —  const.  Inte- 


4  This,  however,  does  not  mean  that  the  temperature  gradient  in  the  CL  is  small 
(see  below). 
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grating  this  equation  over  the  CL  thickness  we  get 


— ■A-eff 


dT 

dx 


+  ^eff 


dT 

dx 


0 


where  (5),  (13)  and  (14)  are  taken  into  account. 

In  the  CL  ere  is  several  orders  of  magnitude  larger,  than  ap. 
Thus,  l/ere  on  the  right  side  of  (15)  can  be  neglected.  Physi¬ 
cally,  Joule  heat  in  the  active  layer  is  generated  mainly  in  the 
electrolyte  phase. 

The  left  side  of  (15)  is  the  sum  of  the  two  one-sided  heat 
fluxes  leaving  the  CL  in  the  opposite  directions:  the  first  one  is 
directed  to  the  backing  layer  and  the  other  to  the  membrane  (Fig. 
1).  Our  goal  is  to  derive  expression  for  the  total  heat  flux  from 
the  CL.  Clearly,  the  total  flux  is  a  sum  of  one-sided  fluxes  and 
we,  therefore,  may  write 


-■A-eff 


dT 

dx 


tot 


(16) 


Eq.  (16)  is  the  desired  boundary  condition  for  external  problems. 
The  sign  “minus”  on  the  left  side  means  that  the  temperature 
drops  along  x  and  the  heat  flux  is  directed  towards  the  backing 
layer  (Fig.  1);  this  sign  should  be  chosen  in  accordance  with  the 
formulation  of  the  external  problem. 

The  first  term  on  the  right  side  of  (16)  represents  the  heat 
flux  due  to  the  electrochemical  reaction.  The  second  term  is  the 
heat  flux  due  to  Joule  heating  in  the  proton-conducting  phase. 
Generally,  all  terms  in  (16)  are  of  the  same  order  of  magnitude 
and  no  one  can  be  neglected. 

The  coefficient  1/3  in  the  second  term  on  the  right  side  of 
(16)  arises  due  to  linear  dependence  jp(x)  (13).  In  the  case  of 
jo/j*  1,  jp(x)  strongly  deviates  from  linear  law.  The  respec¬ 
tive  coefficient  in  Eq.  (16)  can  then  be  calculated  using  the  rela¬ 
tions  [8].  Note  also,  that  if  jo/j*  L  overpotential  is  no  longer 
constant  across  the  CL  and  integration  of  Eq.  (10)  with  rj(x)  from 
Ref.  [8]  leads  to  a  more  complicated  expression.  Physically,  in 
a  CL  with  poor  proton  conductivity  peak  of  reaction  rate  at  the 
membrane  interface  would  lead  to  a  larger  overheat,  especially 
when  the  thermal  conductivity  of  the  CL  is  not  large. 

It  should  be  noted  that  the  Joule  term  in  (16)  is  proportional 
to  the  square  of  local  current  density.  The  distribution  of  lo¬ 
cal  current  in  a  cell  is  usually  very  non-uniform  [14]  and  we 
may  expect  the  effects  due  to  temperature  variation  along  the 
cell  surface.  Furthermore,  in  stacks  with  poor  heat  management 
these  thermal  non-uniformities  may  further  enhance.  Though 
temperature  variation  across  the  CL  in  the  stack  is  still  small, 
the  absolute  temperatures  at  different  points  of  the  cell  surface 
may  differ  quite  strongly.  The  respective  temperature  fields  and 
effects  should  be  studied  with  the  models  of  a  higher  dimen¬ 
sionality.  These  models  may  utilize  the  boundary  condition  (16), 
where  T,  ij  and  jo  should  be  considered  as  local  values. 


4.  Conclusions 


We  consider  the  sources  and  transport  of  heat  in  the  catalyst 
layers  of  low-temperature  polymer  electrolyte  fuel  cells.  Heat 


balance  equation  contains  the  terms,  which  describe  the  ther¬ 
modynamic  heat  generated  or  lost  in  the  electrochemical  reac¬ 
tion,  the  “irreversible”  heat  due  to  transport  of  protons  through 
the  electric  double  layer  at  the  catalyst  particle/electrolyte  inter¬ 
face  and  the  Joule  heat  dissipated  in  the  membrane  phase  due  to 
proton  current. 

The  exact  solution  to  the  heat  transport  equation  shows  that  in 
typical  conditions  the  temperature  variation  across  the  catalyst 
layer  is  negligibly  small.  The  heat  flux  from  the  layer  is,  how¬ 
ever,  not  small.  If  proton  current  density  varies  linearly  with  the 
distance  across  the  layer,  the  heat  balance  equation  can  be  inte¬ 
grated  once  to  yield  the  expression  for  the  heat  flux  generated 
in  the  CL.  This  expression  can  be  used  as  a  boundary  condi¬ 
tion  for  “external”  problems,  which  treat  CL  as  an  infinitely  thin 
interface.  The  limits  of  validity  of  the  resulting  expression  are 
discussed. 


Appendix  A.  Solution  to  heat  transport  equation 


In  this  section  we  return  to  the  general  Eq.  (10).  Above  we 
have  shown  that  heat  in  the  CL  is  generated  mainly  in  the  elec¬ 
trolyte  phase.  Therefore,  the  term  j2 /ere  in  (10)  can  be  neglected. 

Linear  shape  of  proton  current  density  in  the  CL  (13)  is  equiv¬ 
alent  to  the  constancy  of  the  rate  R  across  the  layer  [8];  thus  Eq. 
(5)  reduces  to  jo  =  Rl.  Taking  into  account  this  relation  and  Eq. 
(13)  we  may  rewrite  Eq.  (10)  as 


d 

dx 


jo 

l 


(17) 


Assuming  that  keff  is  constant  and  introducing  dimensionless 
variables: 


*  ~.  jo 

*=7>  Jo  =  — 

'  J* 


Eq.  (17)  transforms  to 


dzT 

dx2 


=  T 


nFrj 

~Xs 


(18) 


A  Sj*l 
nFXe  ff 


70  + 


7o  (1  -  x)2 

(19) 


The  left  side  of  (19)  has  dimension  of  temperature  (K);  there¬ 
fore,  the  expression  in  square  brackets  on  the  right  side  is  dimen¬ 
sionless  parameter: 

jJAS  2apbAS 

a  =  -  =  — - -  (20) 

nFXeff  n  FXef{ 


where  Eq.  (12)  is  used. 

Introducing  the  dimensionless  temperature  and  overvoltage: 


(21) 


where 


7* 


nFb 
A ~S 


Eq.  (19)  transforms  to 


d2f 

dx2 


=  otjo{T  +  rj)  +  2a)o(l  -  xf 


(22) 


1240 


A.A.  Kulikov  sky  /  Journal  of  Power  Sources  162  (2006)  1236-1240 


Table  A.l 

Parameters  for  the  calculations 


Parameter 

Value 

Reference 

Proton  conductivity  of  the  electrolyte  phase, 

i 

[2] 

Op  (Sm  *) 

Tafel  slope,  b  (V) 

0.05 

[2] 

Entropy  change  in  the  ORR,  AS 

326.36 

[13] 

(J  mol'1  KX1) 

Thermal  conductivity  of  water,  A.w 

0.58 

[15] 

(Wn-r1  Kr1) 

Current  density,  jo  (Am  2) 

104 

Assumed 

Overpotential,  r]  (V) 

0.5 

Assumed 

To  qualitatively  illustrate  the  solutions  to  Eq.  (22)  we  will 
assume  that  heat  flux  to  the  membrane  is  negligible,  i.e., 
d T /9x|i=o  =  0;  heat  flux  at  x  =  1  (GDL  interface)  is  then  the 
total  heat  flux  from  the  CL.  At  x  =  1  we  fix  the  temperature: 
T(\)  —  T\ .  Note  that  the  actual  value  of  T\  should  be  determined 
from  solution  of  an  external  problem. 


Solution  to  (22)  subject  to  these  boundary  conditions  is 

qr/—\  (7)  +  fj)COS(wx)  _  2~. 

T(x)  = - ri-  2(1  —  x)  jo 

cos  co 

4  jo[co  sin(&>(l  —  x))  —  cos(&>x)  +  cos  co] 


co~  cos  to 


(23) 


where 


Plot  of  (23)  with  the  parameters  from  Table  A.  1  is  shown  in 
Fig.  A.l.5  Temperature  variation  across  the  CL  is  less  than  0. 1  K 
and  can  safely  be  ignored.6  However,  temperature  gradient  at 
x  —  1  is  not  small.  This  gradient  (given  by  Eq.  (16))  determines 
the  heat  flux  from  the  CL  to  the  backing  layer. 
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5  Parameter  co  is  small;  for  numerical  calculations  it  is  useful  to  expand  (23) 
in  terms  of  co  retaining  only  linear  term. 

6  CCL  generates  water;  thus  we  may  expect  that  in  typical  situations  thermal 
conductivity  of  the  CCL  does  not  differ  much  from  that  of  water  and  for  the  es¬ 
timate  we  put  A.eff  —  Aw.  In  dry  catalyst  layer,  due  to  lower  thermal  conductivity 
the  temperature  variation  can  be  larger. 


